Order and chaos in a triaxial galaxy model with a dark halo component 
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Abstract 

We study the regular or chaotic nature of orbits in a 3D potential describing a triaxial galaxy surrounded 
by a spherical dark halo component. Our numerical calculations show, that the percentage of chaotic orbits 
decreases exponentially, as the mass of the dark halo increases. A linear increase of the percentage of the 
chaotic orbits was observed, as the scale length of the dark halo component increases. In order to distinguish 
between regular and chaotic character of orbits, we use the total angular momentum L tot , as a new indicator. 
Comparison of this new dynamical parameter, with other, previously used chaos indicators, shows that the 
Ltot gives very fast and reliable results in order to detect the character of orbits in galactic potentials. 
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1. Introduction 

We shall study the motion in a 3D composite galactic model described by the potential 

V t {x, y, z) = V g (x, y, z) + V h (x, y, z) , (1) 

where 

v g (x, v,z) = y ln ^ ~ Xx " + ay2 + h7? + c '] ' (2) 

while 

-M h 

Potential (2) describes a triaxial elliptical galaxy with a bulge and a small asymmetry introduced by the 
term —Ax 3 , A << 1 (see Binney & Tremaine, 2008). The parameters a and b describe the geometry of the 
galaxy, while ct is the scale length of the bulge of the galaxy. The parameter v is used for the consistency 
of the galactic units. To this potential, we add a spherically symmetric dark halo, described by the potential 
(3). Here and c/j are the mass and the scale length of the dark halo component respectively. 

The aim of this article is: (i) To investigate the character of orbits in the potential (1) and to determine 
the role played by the dark halo, (ii) To introduce, use and check a new and fast detector, the total angular 
momentum L tot , in order to obtain a reliable criterion to distinguish between order and chaos in galactic 
potentials. 

The Hamiltonian to the potential (1) writes 

H=^( P 2 x + P l+ P 2 z )+V t {x,y,z)=h 3 , (4) 
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Figure 1: (a-d): The (x,p x ) phase plane, when (a, upper left): = 0,/i2 = 516, (b, upper right): = 10000, h,2 = —226, 
(c, lower left): M h = 20000, h 2 = -1007 and (d, lower right): M h = 30000, h 2 = -1788. The values of all other parameters 
are given in text. 



where p Xl p y and p z are the momenta per unit mass conjugate to x, y and z respectively, while ft. 3 is the 
numerical value of the Hamiltonian. 

In this article, we use a system of galactic units, where the unit of length is 1 kpc, the unit of mass 
is 2.325 x 10 7 M Q and the unit of time is 0.97748 x 10 8 yr. The velocity unit is 10 km/s, while G is 
equal to unity. The energy unit (per unit mass) is 100 (km/s) 2 . In the above units we use the values: 
v a = 15, q, = 2.5, a = 1.5, b = 1.8, A = 0.03, while Mh and c/j are treated as parameters. 

2. Results for the 2D system 

The total angular momentum for a star of unit mass, moving in a 3D orbit is 



where L x L y and L z are the three components of the angular momentum along the x, y and z axis given by 




(5) 



L 



X 



= yz-yz 



f ~ J qj Z^C ZOO j 

L z = xy~ xy, 



(6) 



where the dot indicates derivative with respect to the time. For a 2D system we set in (6), z = z = 0, that is 
L to t reduces to L z . In this research, we shall use the plot of the L to t vs time in order to distinguish regular 
from chaotic motion. 




Figure 2: (a-d): The (x,p x ) phase plane, when = 10000 and (a, upper left): = 10.5, /12 = —135, (b, upper right): 
Ch = 13, /12 = —55, (c, lower left): = 15.5, /12 = H and (d, lower right): = 18, /12 = 68. The values of all other 
parameters are given in text. 




Figure 3: (a-b): (a, left): A plot of the area A% covered by chaotic orbits versus and (b, right): A plot of the area A% 
covered by chaotic orbits versus c^. The values of all the other parameters are given in text. 



Our next step is to study the properties of the 2D dynamical system, which comes from potential, (1) if 
we set z — 0. The corresponding 2D Hamiltonian writes 

H 2 = 1 -(pl+p 2 y )+V t (x,y) = h 2 , (7) 

where h 2 is the numerical value of the Hamiltonian. We do this in order to use the results obtained for the 
2D model in the study of the more complicated 3D model, which will be presented in the next Section. 
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Figure 4: (a-b): (a, left): A plot of the L.C.E versus and (b, right): A plot of the L.C.E versus c^. The values of all the 
other parameters are given in text. 




Figure 5: (a-d): (a, upper left): An orbit in the 2D potential, (b, upper right): The corresponding L.C.E, (c, lower left): The 
P(f) indicator and (d, lower right): The Ltot indicator. The motion in chaotic. See text for details. 

Figure la-d shows the (x,p x ), (y — 0,p y > 0) phase plane, for four different values of the mass of the 
dark halo. The values of all the other parameters are: vo = 15, q, = 2.5, a = 1.5,6 = 1.8, A = 0.03 and 
Ch = 8. Fig. la shows the phase plane, when the system has no halo component, that is when Mh = 0. The 
value of h 2 is 516. Fig. lb is similar to Fig. la but when Mh = 10000 and h 2 — —226. Fig. lc is similar 
to Fig. la but when M h = 20000 and h 2 = -1007. Fig. Id shows the phase plane when M h = 30000 and 
h 2 = -1788. 

Fig. 2a-d is similar to Fig. la-d but when Mh is 10000, while Ch is treated as a parameter. All the other 
parameters are as in Fig. 1. In Fig. 2a we have Ch = 10.5 and h 2 = —135. In Figure 2b the values of Ch and 
h 2 arc 13 and -55 respectively. In the phase plane shown in Fig. 2c, we have taken Ch = 15.5 and h 2 = 11. 
In Fig. 2d we have chosen Ch — 18 and h 2 = 68. 

Figure 3a shows the percentage of the phase plane A% covered by chaotic orbits as a function of the 



4 




mass of the dark halo, for two different values of x max . The values of the parameters are: v = 15, Cf, = 
2.5, a = 1.5, b — 1.8, A = 0.03 and Ch = 8. Figure 3b shows a plot between A% and Ch- The values of the 
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Figure 8: (a-d): Similar as Fig. 5a-d but for a chaotic orbit. 
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Figure 9: (a-b): Similar to Fig. 4a-b but for the 3D potential. The values of all other parameters are given in text. 

parameters are: v = 15, c& = 2.5, a = 1.5, b = 1.8, A = 0.03, Mh = 10000. Figures 4a and 4b show a plot of 
the Lyapunov Characteristic Exponent (L.C.E) versus M h or versus respectively. 

In what follows, we shall investigate the regular or chaotic character of orbits in the 2D Hamiltonian 
(7) using the new dynamical detector L to t ■ In order to see the effectiveness of the new method, we shall 
compare the results with two other indicators, the classical method of the L.C.E and the P(f) spectral 
method, used by Karanis & Vozikis (2007). This method uses the Fast Fourier Transform (F.F.T) of a series 
of time intervals, each one representing the time that elapsed between two successive points on the Poincare 
(x,p x ) phase plane for 2D systems, while for 3D systems they take two successive points on the plane z = 0. 

Figure 5a shows an orbit with initial conditions: xq — — 1.0, ya = p x o = 0, while the value of p y o is always 
found from the energy integral for all orbits. The values of all the other parameters and the energy are as in 
Fig. la. One observes in Figure 5b that the L.C.E, which was computed for a period of 10 5 time units, has 
a value of about 0.18 indicating chaotic motion. The same result is shown by the P(f) indicator, which is 
given in Figure 5c. Figure 5d shows a plot of the L tot vs time, for a time interval of 100 time units. We see 
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Figure 10: (a-d): (a, upper left): An orbit in the 3D potential, (b, upper right): The corresponding L.C.E, (c, lower left): The 
P(f) indicator and (d, lower right): The Ltot indicator. The motion in chaotic. See text for details. 

that the diagram is highly asymmetric, with a number of small and large peaks. The above characteristics 
suggest that the corresponding orbit is chaotic. 

Figure 6a shows an orbit with initial conditions: x = 8.8, yo = Pxo = 0. The values of all the other 
parameters and the energy are as in Fig. lb. As we see, this is a quasi periodic orbit. Therefore the 
L.C.E of this orbit goes to zero, as it is clearly seen in Figure 6b. The P(f) indicator in Figure 6c shows a 
small number of peaks, also indicating regular motion. The plot of the L tot given in Figure 6d is now quasi 
periodic, with symmetric peaks, indicating regular motion. 

Figure 7a-d is similar to Fig. 6a-d for an orbit with initial conditions: x = — 9.36,yo = Pxo — 0, while 
the values of all the other parameters and the energy are as in Fig. 2a. As we see, the orbit is quasi periodic 
and this fact is indicated by all three dynamical parameters. On the contrary, the orbit shown in Figure 8a 
has initial conditions: xq = 10, yo = p x o = and the values of all the other parameters and the energy are 
as in Fig. 2d. The orbit looks chaotic and this is indicated by the L.C.E, the P(f) and the L tot shown in 
Figs. 8b, 8c and 8d respectively. 

A large number of orbits in the 2D system were calculated for different values of the parameters. All 
numerical results suggested that the L tot is a fast and reliable dynamical parameter and can be safely used 
in order to distinguish ordered from chaotic motion. 



7 



t do 4 ) 





100 



Figure 11: (a-d): Similar as Fig. lOa-d. The motion is regular. 



3. Results for the 3D system 

The regular or chaotic nature of the 3D orbits is found as follows: we choose initial conditions (xo, p x o, z o), 
Vo = Pzo = 0, such as (x , p x o) is a point on the phase plane of the 2D system. The point (xn,f>xo) lies inside 
the limiting curve 



-pl + V t {x) = h 2 , 



(8) 



which is the curve containing all the invariant curves of the 2D system. We choose h 3 = h 2 and the value 
of p y o for all orbits is obtained from the energy integral (4) . Our numerical experiments show that orbits 
with initial conditions (x a ,p x0 , z a ), y — p za = 0, such as (x a ,p x0 ) is a point in the chaotic regions of Figs, 
la-d and 2a-d for all permissible values of z , produce chaotic orbits. 

Our next step is to study the character of orbits with initial conditions {xo,p x o, zq), yo = Pzo = 0, such 
as (xo,p x0 ) is a point in the regular regions of Figs, la-d and 2a-d. It was found, that in all cases the regular 
or chaotic character of the above 3D orbits depends strongly on the initial value z . Orbits with small of 
zq are regular, while for large values of zq they change their character and become chaotic. The general 
conclusion, which is based on the results derived from a large number of orbits, is that orbits with values of 
z > 0.75 are chaotic, while orbits with values of z < 0.75 are regular. 

Figure 9a shows the L.C.E of the 3D system, as a function of the mass of halo, for a large number of 
chaotic orbits, when Ch = 8. Figure 9b shows he L.C.E as a function of Ch when Mh = 10000. 

Figure lOa-d is similar to Fig. 8a-d but for a 3D orbit. The orbit shown in Fig. 10a looks chaotic. The 
initial conditions are: x = 2.0, p x0 = 0, z = 0.5. Remember that all orbits have y = p z0 = 0, while the 
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value of p y o is always found from the energy integral. The values of all the other parameters and the energy 
/13 are as in Fig. 2b. The L.C.E shown in Figure 10b assures the chaotic character of the orbit. The P(f) 
given in Figure 10c also suggests chaotic motion. The same conclusion comes from the L tot , which is shown 
in Figure lOd. Figure lla-d is similar to Fig. lOa-d but for a quasi periodic 3D orbit. The initial conditions 
are: xo = 5.0, p x o = 0, zq = 0.1. The values of all the other parameters and the energy /13 are as in Fig. lc. 
Here one observes, that all the three detectors support the regular character of orbit. 

The conclusion for the study of the 3D model, is that the L to t detector can give reliable and very fast 
results for the character of the orbits. There is no doubt, that the L tot is faster than the two other indicators, 
used in this research. Therefore, one can say that this indicator is a very useful tool for a quick study of the 
character of orbits in galactic potentials. 

4. Discussion and conclusions 

The main conclusions of this research arc the following: 

1. The percentage of the chaotic orbits decreases as the mass of the spherical halo increases. Therefore, 
the mass of the dark halo can be considered as an important physical quantity acting as a controller of 
chaos in galaxies showing small asymmetries. 

2. One expects to observe a smaller fraction of chaotic orbits in asymmetric triaxial galaxies with a 
dense spherical halo, while the fraction of chaotic orbits would increases in asymmetric triaxial galaxies 
surrounded by less dense spherical dark halo components. 

3. It was found that the L.C.E in both the 2D and the 3D models decreases, as the mass of halo 
increases, while the L.C.E increases as the scale length c h of the halo increases. This means that not only 
the percentage of chaotic orbits, but also the degree of chaos is affected by the mass or the scale length of 
the spherically symmetric dark halo component. 

4. The Ltot gives fast and reliable results regarding the nature of motion, both in 2D and 3D galactic 
potentials. For all calculated orbits the results given by the L to t coincide with the outcomes obtained using 
the L.C.E or the P(f) spectral method. The advantage of the L tot is that is faster than the above two 
methods. 
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